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Abstract 

Monte Carlo (MC) simulations and series expansions (SE) data for the 
energy, specific heat, magnetization, and susceptibility of the three-state and 
four-state Potts model and Baxter- Wu model on the square lattice are an- 
alyzed in the vicinity of the critical point in order to estimate universal 
combinations of critical amplitudes. We also form effective ratios of the 
observables close to the critical point and analyze how they approach the 
universal critical- amplitude ratios. In particular, using the duality relation, 
we show analytically that for the Potts model with a number of states < 4, 
the effective ratio of the energy critical amplitudes always approaches unity 
linearly with respect to the reduced temperature. This fact leads to the 
prediction of relations among the amplitudes of correction-to-scaling terms 
of the specific heat in the low- and high-temperature phases. It is a com- 
mon belief that the four-state Potts and the Baxter- Wu model belong to 
the same universality class. At the same time, the critical behavior of the 
four-state Potts model is modified by logarithmic corrections while that of 
the Baxter- Wu model is not. Numerical analysis shows that critical ampli- 
tude ratios are very close for both models and, therefore, gives support to 
the hypothesis that the critical behavior of both systems is described by the 
same renormalization group fixed point. 
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The fixed points of tlie renormalization group define the universal behav- 
ior of a system through a set of critical exponents and universal combinations 
of critical amplitudes The universality concept divides all systems at crit- 
icality into a number of universality classes. It is instructive to know the set 
of values of the critical exponents and of the universal combinations of critical 
amplitudes for a given universality class. 

The two-dimensional Potts model 0] is the simplest model which exhibits 
a phase transition. It is solved exactly at the critical point for any number of 
spin components q and it is known that for g < 4 it has a continuous phase 
transition while for g > 4 the phase transition is of the first order. The model 
has a great theoretical interest as new theories may be tested in this model. 

At the same time, these models may have some practical interest as they 
may be realized in an adsorbate lattice placed onto a clean crystalline sur- 
face. The full classification of such systems with continuous transitions is 
known theoretically [3|. There are experiments in which some of them real- 
ize the 3-state and 4-state Potts models [3] and the critical exponents can be 
experimentally estimated. 

Critical exponents for the Potts model with g < 4 can be computed 
exactly by different theoretical techniques j^, 0] . The values of the thermal 
critical exponents and of the magnetic critical exponents follow from the 
identification of the dimensions of the conformal algebra operators 

Nowadays, there is no doubt on the values of the leading critical expo- 
nents whereas the values of the correction-to-scaling exponents are still under 
discussion, as well as the values of the universal ratios of the critical ampli- 
tudes. Our presentation is intended to give a short review of the research on 
the subject. 

The Potts model Hamiltonian (see review [?!] for details) can written as 
H = — ^SiSj , where is a spin variable taking integer values between 
and q—l, and the sum is restricted to the nearest neighbor sites {ij) on the 
square lattice. 

Close to the critical temperature Tc at which the continuous phase tran- 
sition occurs, the residual magnetization M and the singular part of the 
reduced susceptibility x and of the specific heat C of the system in zero ex- 
ternal field are characterized by the critical exponents /?, 7, and a and by 
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the critical amplitudes B, r±, and A± 



M{t) ^ B{-Tf, r < (1) 
X±(r) ^ r±|r|-^ (2) 

C±{r) ^ ^|r|-, (3) 
a 

where r is the reduced temperature r = (T — Tc)/T and the labels ± refer to 
the high-temperature and low-temperature sides of the critical temperature 
Tc. The critical amplitudes are not universal by themselves but some com- 
binations of them, f.e., Aj^/A^, r_|_/r_, and T^A^/ B"^, are universal pj] due 
to the scaling laws. 

On the square lattice, in zero field, the model is self-dual. The duality 
relation 

(e/3 - 1) (e/^* -l)=q (4) 

fixes the inverse critical temperature to /5c= In(H-yg). The values E{P) and 
E{(3*) of the internal energy at dual temperatures are simply related through 

(l-e-^)E(/3) + (l-e-^*)E(/5*) = 2. (5) 

Dual reduced temperatures r and r* can be defined by /3=/3c(l— t) and 
/3*=/3c(l+T*). Close to the critical point, r and r* coincide through linear 
order, since T*=T+'^^^^^±f^THO{T^). 

The ratio of the free energy critical amplitudes Aj^/A_ is equal to unity 
due to duality. Moreover, duality relations may be used to understand the 
dependence on temperature of the effective amplitude functions which may 
be constructed from the energy in the symmetric phase E^{t) and in the 
ordered phase E^{t*) 

A+(r) = a{l-a)l3,{E+{T)-E^)T'^-\ (6) 
A_(r*) = a{l-a)(3,{Eo~E_{r*)){T*r-' (7) 



as an effective amplitude ratio 

A+(r) {E^{t) - E,)t-'^ 



a— 1 ' 



(8) 



where the constant Eq is the value of the energy at the transition temperature, 
E^=E{[3,)=-l-l/^. 
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Table 1: Exact values of critical exponents and ratios of critical amplitudes for the Ising 
model (2-state Potts model). 





a 


/5 


7 




r+/r^ 


= T+A+/B'' 


1 





1/8 


7/4 


1 


37.69365... 


0.318569... 



Evaluating expression (IHl) for small r and denoting aq = —EoPc^ 
~\y^\ we obtain 



ln(l+^) 

, VV^ W U U CLIXJ 

A, 



AJt*] 



l + (3-a)a,r + 0(r^+"). (9) 



Note the linear dependence on r of the effective amplitude ratio. 

The 2-state Potts model is equivalent to the Ising model which was solved 
exactly Q] (see Ref. ^ for details). The susceptibility behavior was under- 
stood in the paper by Wu, McCoy, Tracy and Barouch [l(|. It turns out 
that there exist only integer corrections to scaling (for a recent and detailed 
discussions we refer readers to Refs. 11, 12, 3]). Values of the critical ex- 



ponents and of some amplitude ratios are presented in Table [TJ 
The critical behavior of the susceptibility reads as 

X{r) = r±|r|-^^,,„,(|r|^) + y^ir), (10) 

where A'corr (|t|^) is the correction-to-scaling function and 3^w(t) represents 
an analytic expression ( "background term" ) which accounts for non-singular 
contributions to susceptibility. 

Set of values of the thermal critical exponents x^^ and of the magnetic 
critical exponents Xa„ are known analytically ^ ^ 

_n^_+ny _ {2n - if - 

2-2/ ' 4(2 -y) ^ ^ 

in terms of the variable y linked to the number of states q by cos ^ = \\[q- 
For the 3-state Potts model there is a finite number of correction terms [6| , 
= 14/5, = 6, and x„^ = 4/3. The leading correction-to-scaling 
contribution is A = — (2 — x^^)/ {2 — x^J = 2/3 and it was first supported in 



numerical simulation 14 



Clear evidence for these leading correction to the scaling behavior may 
be seen in Figure [Tj where we plot the difference of the longitudinal and 
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Table 2: Exact values of critical exponents and estimates of the ratios of critical amplitudes 
for the 3-state Potts model. 



V 


a 


/? 


7 


r+/ri 


Tt/Fl 


i?J=r+A+/i?^ 


Remark 


5/6 


1/3 


1/9 


13/9 








Exact result 










13.848 


0.327 


0.1041 


[15, 16] 










13.83(9) 


0.325(2) 




[19] - SE 










13.83(9) 


0.3272(7) 


0.1044(8) 


[18] - SE 










13.86(12) 


0.322(3) 


0.1049(29) 


[18] - MC 



transverse susceptibilities (xl — Xr) (note the cancelation of background 
terms fIlOp in the difference) multiplied by the leading behavior factor |rp, 
as a function of |Tp/^. The ratio of the susceptibilities {xt/xl) is shown 
in Figure [2] and may be used to estimate the universal ratio of associated 
amplitudes r^/r^. 

Analytical predictions for the amplitude ratios in Potts model for q = 
2,3, and 4 were given in the papers |l5l . Il6l |. The values are shown in the 



Table [2] together with numerical estimations from Monte Carlo (MC) and 



series expansions (SE) analyses [l7|, [l8|, [l9 



The coincidence of the data is 

to 
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a good indication for the validity of both two-kink approximation [15| ^ 

the exact scattering theory i^j for 3-state Potts model and of the analysis 
of MC data and series expansions. 

The analysis of the 4-state Potts model is much more complicated be- 
cause in addition to the corrections to scaling there are confluent logarithmic 



corrections |2l|, |22j. The result of the analysis of MC data and SE 23l. 24l . |19 



is shown in the Table [3] together with the analytical estimates 15|, Il6 |. 

Conformal field theory predicts the set of renormalization group (RG) 
exponents y^^=2—x^^=3/2, 0, -5/2, ... . The leading correction-to-scaling 
exponent y^^ vanishes and gives rise to a logarithmic behavior 2l|. In our 
recent publication [il, 2^1 we revised the renormalization group equations 
and included in our analysis the known form of the logarithmic corrections 
and of the next-to-leading corrections, taking into account the width of the 
temperature region window examined. The set of magnetic exponents for 
the Potts model Xo-„=l/8, 9/8, 25/8, ... translates into the magnetization 
exponent (3=1/12 and the leading correction-to-scaling exponent Ao-=2/3. 
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Finally, the following behavior of the susceptibility is assumed 



X+{r) =T+r-''/'^g^/\-\nT){l + a+T^/^ + b+T) + D+, (12) 

where the function Q contains a universal correction function £ 22, 24, 2^ 
and the leading nonuniversal correction function JF 



g{-\n\T\) 
^(-ln|r|) 



ln|r|) X S{-\n\T\) x J^(-ln|r| 
31n(-ln|r|)~ 
4' 



1 + - 



In r 



^ 31n(-ln|r|) 
4 -Inlrl 



l + - 



^(-ln|r|) 



1 



Ci C2ln(-ln|r|) 



In |r| 



ln|r|)2 



4(-ln|r|) 
1 



(13) 

(14) 
(15) 



We fit our data to estimate the amplitude r+, the coefficient of the leading 
correction to scaling a+ in Eq. f|T2l) . and coefficients Ci and C2 in Eq. f|T5|) . 

It is obvious that the logarithmic corrections (the whole function Q{— In |r|)) 
cancels in simple ratios, like A+/A_, T+/Tl, Tt/Tl, etc. This has been 
demonstrated analytically for the effective ratio A+/A_ (see Eq. [2]). We note 
also that the RG analysis predicts powers of logarithmic corrections 

to specific heat a'=—l, susceptibility 7'=3/4, and magnetization /?'=— 1/8 
such that they cancel in all universal ratios. For example, the universal 
amplitude ratio may be calculated as the limit of the ratio of functions 



R. 



c 



(E_(M)-i^o)x-(M) . 

r^o M(|r|)2 ^ ^^'^ 



(16) 



where Eq = E{0) = V2. One can check that in the ratio (|T6l) not only 
powers of |r| cancel but also powers of S. In the ratio, the magnetization M 
and the energy difference E'-drl) — Eq have only singular contributions and 
the only systematic deviation may come from the background correction to 
susceptibility x_(|r|). It was shown in ISj that the contribution from this 
background correction is negligible in the critical region window, and the 
estimator ( |T6l) tends to the value 0.0055(1) as r — » 0. 

The Baxter- Wu [26] model is defined on a triangular lattice, with spins 
CTj located at vertices. The three spins forming a triangular face are coupled 
with a strength J, and the Hamiltonian reads 
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Figure 1: 3-statc Potts model. Difference of susceptibilities {xl — Xt)\t\'^ as function of 
|^|2/3_ rpjjg almost linear dependence supports the value 2/3 for the power of the leading 
correction to scaling. 
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Table 3: Exact values of critical exponents and estimations of the ratios of critical ampli- 
tudes for the 4-state Potts model. 



V 


a 


/? 


7 


r+/ri 


Tt/Fl 




Remark 


2/3 


2/3 


1/12 


7/6 








Exact result 










4.013 


0.129 


0.00508 


[15, 16] 










3.5(4) 


0.11(4) 




[19] - SE 










3.14(70) 




0.0068(9) 


[23] - MC 










6.93(6) 


0.1674(30) 


0.00512(13) 


[24, 25] - MC 










6.30(1) 


0.1511(24) 


0.00531(5) 


[24, 25J - SE 



= - J ^ (Tiajak, (17) 

faces 

where the summation extends over all triangular faces of the lattice, both 
pointing up and down. The ground state is four-fold degenerate and the 
critical exponents are found to be the same as for the 4-state Potts model. 
The exact behavior of the magnetization, energy, and specific heat are 



known analytically |26l. l27l . 1281 ]. An analysis of Monte Carlo data was per- 
formed by two of us [29| and preliminary estimates shows that the values of 
the susceptibility amplitude-ratio T^/Fx, ^ 6.9 and of the ratio ~ 0.005 
are very near to those obtained from our analysis of MC and SE data for 
the 4-state Potts model (see Table [3]). We have to note that logarithmic 
corrections to scaling are absent in the critical behavior of Baxter- Wu model 
and this gives us more confidence in our analysis. 



Delfino and Grinza [30j use the same analytical approach as in [15| to 
study the Ashkin- Teller model which also belongs to the 4-state Potts model 
universality class with some particular choice of coupling constants. This 
leads to the estimatation r+/rL ~ 4.02 and Pt/Pl ~ 0.129. This result is 
also very near to those for 4-state Potts model (see second entry in Table [31 ) 

A possible explanation of the deviation of our results from the analytical 
predictions may be explained as follows: the two-kink approximation is exact 
for the 2-state Potts model, it gives good accuracy for the 3-state Potts model, 
but it may be insufficient to produce accurate values for 4-state Potts model. 
Further analyses have to be done to resolve the contradiction among these 
results. 
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Figure 2: 3-state 
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